EUROPEAN LABORATORY FOR PARTICLE PHYSICS 



CERN-EP/99-045 
March 19, 1999 



Test of the Flavour Independence of 
as using Next-to-Leading Order 
Calculations for Heavy Quarks 



The OPAL Collaboration 



Abstract 

We present a test of the flavour independence of the strong couphng constant for 
charm and bottom quarks with respect to hght (uds) quarks, based on a hadronic 
event sample obtained with the OPAL detector at LEP. Five observables related to 
global event shapes were used to measure Ug in three flavour tagged samples (uds, 
c and b). The event shape distributions were fitted by 0{a1) calculations of jet 
production taking into account mass effects for the c and b quarks. We find: 

a^/a^ds ^ 0.997 ± 0.038 (stat.) ± 0.030 (syst.) ± 0.012 (theory) and 
ab/Q,^ds ^ 0.993 ± 0.008 (stat.) ± 0.006 (syst.) ± 0.011 [theory). 
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1 Introduction 



According to the theory of strong interactions, Quantum Chromodynamics (QCD) |Jl|], 
the strong couphng constant is the same for all quark flavours. Therefore a precise 
measurement of as for the individual quark flavours is an important test of this theory. 
The coupling constant for charm and bottom quarks can be compared to for light 
(uds) quarks by measuring the ratios a^/a"'^^ and a^/a"'^^ using hadronic events of the 
type e"'"e~ cf^g. The coupling constants are measured using an event sample originating 
from a pair of light quarks (uu, ss or dd), c quarks (cc) or b quarks (bb), respectively. If 
either of the above ratios deviates significantly from unity then this may indicate physics 
beyond the Standard Model. 

In many QCD studies at LEP corrections due to quark mass effects can be safely 
ignored because they typically appear as powers of the ratio of the quark mass to the 
total energy. However, in those studies where is determined in event samples enriched 
in heavy quarks, mass effects become non- negligible. Gluon emission by bottom quarks, 
and to a lesser extent charm quarks, will be suppressed largely due to the reduced phase- 
space available. Observables sensitive to the three-jet rate measured in heavy quark events 
will be modified with respect to the same quantities measured in light quark events. 

Tests of the flavour independence of as have previously been conducted at both LEP 
and SLC 0-0] using ratios of as for one flavour over as of either the complementary [] 
or inclusive quark mixture. In all these tests flavour independence of was confirmed. 
These tests, however, used either massless QCD calculations or a leading order calculation 
with massive quarks The leading order calculation includes the process e"'"e~ qqS'fi', 
with mass effects, but virtual corrections to the process e"'"e~ — > qq(? are not included. 

Complete next-to-leading-order (NLO) calculations of the heavy-flavour production 
cross section in e"'"e~ collisions, including quark mass effects, have been published re- 
cently p-pT]]. A comparison between these calculations has been made and they were 



found to be in agreement |TT|. In recent publications DELPHI [jT2[ and SLD [0 re- 
port on measurements of al/a'f'^ where NLO massive calculations ||, were used to 
account for mass effects in b quark events. These results are consistent with the flavour 
independence of as- 

In this study we use the results of P. Nason and C. Oleari along with theoretical 
predictions assuming massless quarks in fits to global event shape distributions in order to 
determine a^^^ , a^^/a^"^^ and a^/a^^^. The results presented here are intended to update 
and supersede the corresponding OPAL results in [Q, Q insofar as we now have greatly 
increased the charm event statistics and have used improved theoretical predictions. 

This paper is organised as follows. In Sect. ^ the parts of the OPAL detector most 
important to this analysis are described. In Sect. |^ the hadronic event sample and the 
Monte Carlo event sample are introduced. In Sect. ^ and ^ the flavour tagging methods 
are described. In Sect. |] the event shape observables used in this study are introduced 
and the procedure for correcting the event shape distributions is explained. Next, in 
Sect. 0, the procedure for fitting the NLO QCD prediction to the corrected distributions 
is explained and in Sect. IS the systematic uncertainties that have been taken into account 



^For example, as for c quarks was compared to as for a mixtm'e of u, d, s and b quarks. 
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are discussed. In Sect. ^, the results of the test of flavour independence of are presented. 
Finally, in Sect. |l^, conclusions are drawn. 



2 The OPAL detector 

The OPAL detector operates at the LEP e^e~ collider at CERN. A detailed description 



can be found in Refs. [H, 15|. The analysis presented here relies mainly on the recon- 
struction of charged particle trajectories and momenta in the central tracking chambers, 
on energy deposits ("clusters") in the electromagnetic calorimeters and on information 
from the silicon micro-vertex detector. 

All tracking systems are located inside a solenoidal magnet which provides a uniform 
magnetic field of 0.435 T along the beam axi^. The magnet is surrounded by a lead glass 
electromagnetic calorimeter and a hadron calorimeter of the sampling type. Outside the 
hadron calorimeter, the detector is surrounded by a system of muon chambers. There are 
similar layers of detectors in the barrel (|cos^| < 0.82) and endcap (|cos^| > 0.81) regions. 

The central tracking detector consists of a silicon micro- vertex detector [|1^] and three 
drift chamber devices: the vertex detector, a large jet chamber, and surrounding z- 
chambers. The silicon micro-vertex detector, close to the beam pipe, consists of two 
layers of silicon strips with a single-hit resolution of about 7/im in the r0 plane. The 
vertex chamber is a cylindrical drift chamber covering a range of | cos^| < 0.95. Its single 
hit resolution is 50 /im in the r0 plane and 700 fim in the z direction. The jet chamber is 
a cylindrical drift chamber with an inner radius of 25 cm, an outer radius of 185 cm, and a 
length of about 4 m. Its spatial resolution is about 135 /im in the r0 plane from drift time 
information and about 6 cm in the z direction from charge division. The z-chambers pro- 
vide a more accurate z measurement with a resolution of about 300 fim. In combination. 



the three drift chambers yield a momentum resolution of (Jpt/pt ~ y 0.02^ -|- (0.0015 ■ pt)"^ 
for I cos(^)| < 0.7, where pt is the transverse momentum in GeV/c. 

Electromagnetic energy is measured by lead glass calorimeters surrounding the solenoid 
magnet coil. They consist of a barrel and two endcap arrays with a total of 11704 lead 
glass blocks covering a range of | cos6\ < 0.98. 



3 Event sample and Monte Carlo simulation 

This analysis is based on a sample of 4.4 million hadronic decays of the Z° recorded 
with the OPAL detector between 1990 and 1995. Hadronic Z° decays were selected by 
placing requirements on the number of reconstructed tracks and the energy deposited 



in the calorimeter. A detailed description of the criteria is given in |T6|. The parts 
of the detector essential for the present analysis (central detector and electromagnetic 
calorimeter) were required to be fully operational. The track selection criteria were the 



same as presented in a previous OPAL study |T^. The number of accepted tracks was 



^In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis 
points upwards and the z axis points in the direction of the electron beam. The polar angle 9 and the 
azimuthal angle are defined with respect to the z- and a;-axes, respectively, while r is the distance from 
the z-axis. 
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required to be at least five to reduce r+r^ background. Clusters of electromagnetic energy 
were used if their observed energy was greater than 0.25 GeV, and known noisy channels 
in the detector were ignored. The event thrust axis was determined using all accepted 
tracks and clusters, and its direction was required to fulfil the condition | cos6'Th| < 0.9 in 
order that the event be well contained. Using these selection criteria, Monte Carlo studies 
indicate that, within the chosen range of cos6'Th; 99.86±0.07% of hadronic Z" decays are 
accepted, with a contamination of about 0.14% from t~^t~ events, and around 0.07% from 



two-photon interactions [17 



To correct the measured event shape distribitions (see Sect. |6.1| ), 4 million hadronic 
decays of the have been simulated using the JETSET 7.4 Monte Carlo model with 
parameters tuned to represent OPAL data well . For all simulated events heavy quark 
fragmentation has been implemented using the model of Peterson et al. [^. All events 
have been passed through a detailed simulation of the OPAL detector [JT| before being 
analysed using the same programs as for data. 



4 Selection of uds and b quark events 

Two event samples, one enriched in uds quark events and one enriched in b quark events, 
were selected by first determining the number of tracks in each event (Nsig) with a large 
impact parameter significance, b/cxh > 2.5. Here b is the distance of closest approach in 
the xy plane of the track to the e~*'e~ interaction point (IP) and at its error. The sign of b 
was determined with respect to the crossing point between the track and the jet axis, b is 
positive if the track crosses the jet axis downstream of the IP and negative otherwise. An 
event was classified as having a high probability to have come from light quarks if Ngig 
was zero and from b quarks if Ngig > 5. In the text below these will be referred to as the 
"uds-tag" and the "b-tag" event samples. 

For the uds and b quark event selection all data recorded during 1994 were used. This 
represents 1.4 million events. The analysis was restricted to only the 1994 data because of 
the uniform configuration of the silicon micro- vertex detector during this time period. In 
addition to the track selection criteria outlined above, the track was required to contain 
at least one silicon hit and an algorithm was applied to reject tracks which were consistent 
with arising from photon conversions In order for the event to be contained within 



the acceptance of the silicon micro- vertex detector the event's thrust axis was further 
restricted to lie within the range | cos6'Th| < 0.7. 

The distribution of Ngig is shown in Fig. |l| compared with the result of Monte Carlo 
simulation. The simulation is decomposed into the contributions from uds, c and b quark 
events. The general agreement between the sum of the three contributions from Monte 
Carlo and the data is good. 

This tagging procedure resulted in 325 111 events selected for the uds-tag event sample 
and 71 521 for the b-tag event sample. The fiavour compositions of the uds-tag and the 
b-tag event samples as determined from the Monte Carlo, along with their combined 
statistical and systematic errors, are presented in Table |I[ The efficiency for tagging uds 
events was about 35% and for tagging b events 23%. 
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4.1 Uncertainties in Flavour Composition 

The systematic errors on the flavour compositions resuh from uncertainties in the detector 
modelhng and imprecise knowledge of physics processes. For each source of uncertainty 
some aspect of the Monte Carlo was varied and the flavour composition of the uds-tag 
and b-tag event samples were recalculated. The difference between the flavour fractions 
calculated for each variation and the central value was taken as a systematic error on the 
determined flavour fractions. 



4.1.1 Detector Modelling Uncertainties 

Since the determination of the flavour fractions of the uds-tag and the b-tag event sample 
depend on Monte Carlo this requires an accurate simulation of the detector resolution for 
charged tracks measured with the silicon micro-vertex detector. The simulation has been 
tuned to reproduce the tracking resolutions seen in data by studying the impact parameter 
distributions of tracks, as functions of track momentum, polar angle and the different sub- 
detectors contributing hits. This tuning procedure was affected by uncertainties in the 
radial alignment within the silicon micro-vertex detector, the efficiency for associating 
silicon hits to tracks and the modelling of known inefficient regions, and the overall track 
reconstruction efficiency. These uncertainties were evaluated according to the procedure 
given in and their effect on the flavour fractions determined from Monte Carlo was 
calculated. In addition, the agreement between data and Monte Carlo in Fig. Q was 
improved by degrading the resolution of impact parameters in the Monte Carlo simulation 
by 5%. This was done by applying a single multiplicative factor j3 to the difference between 
the reconstructed and true impact parameters. Fig. ^ is shown with this smearing applied. 
To evaluate the sensitivity of the determined flavour fractions to the tracking resolution 
the Monte Carlo smearing was removed. 



4.1.2 Physics Modelling Uncertainties 

The tagging efficiencies of the uds and b tags for the various quark flavours are also 
sensitive to various physics input parameters in the Monte Carlo simulation. The rate of 
gluon splitting to cc was varied in the range (2.38 ± 0.48) x 10~^, based on the OPAL 
measurement |2^, and the rate of gluon splitting to bb was varied in the range (3.1 ± 



1.1) X 10~^, based on theoretical expectation [^. The error due to these two sources 
was negligible. The production fractions of the different weakly decaying b hadrons were 
varied according to the experimental uncertainties [^. The production fractions of the 
weakly decaying c hadrons, the fragmentation of b and c quarks, and the charged decay 
multiplicities and lifetimes of b and c hadrons were varied according to the prescription 
given in [Q. The largest effect on the tag flavour fractions comes from the b hadron 
charged decay multiplicity. All these uncertainties, along with the uncertainties due to 
the detector modelling, were added in quadrature and are given as the errors on the 
flavour fraction shown in Table |l|. 
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5 Selection of c quark events 



Events having a high probabihty to have originated from cc were identified by the presence 
of a highly energetic D*"*" mesonQ. These events will be referred to in the text below as the 
"c-tag" event sample. For the D*+ reconstruction and the determination of the fiavour 
composition of the c-tag event sample the methods of a previous OPAL study have been 
used p^- brief, five D*+ decay modes were searched for: 



D07r+ 








U K- 




h 


"3 prong" 


U K- 






"satellite" 


U K- 






"5 prong" 


U K- 


-e+ 




"electron" 


U K- 






"muon" . 



No attempt was made to reconstruct the 7r° in the satellite channel, nor the neutrino 
direction or energy in the electron and muon channels. The last two channels are referred 
to as "semileptonic" channels in the following text. A number of tracks appropriate for 
the selected channel were combined to form a candidate and their invariant mass was 
calculated. Candidates were selected if the reconstructed mass lay within the expected 
range for that channel. After adding a further track as a possible pion from the D*+ decay, 
the combined mass was calculated and the candidate was selected if the mass difference 
AM = M£,*+ — Mdo was within certain limits. Some of the D*"*" selection criteria are 
given in Table |^. 

To reduce the background in the c-tag event sample we required > 0.4, where 
a;D*+ = E^i+ / -Ebeam is the scaled energy^ of the D*"*" meason. For the 5-prong event 
selection we required xd*+ > 0.5. This cut is effective in rejecting D*"*" mesons which 
originate from cascade decays of B hadrons and from events where a gluon splits into a 
pair of charm quarks. 

The distribution of AM for all five decay channels is presented in Fig. The points 
with error bars are the signal candidates and the solid histograms are background estima- 
tor distributions constructed from data. The background estimator was constructed by 
choosing the candidate for the pion in the D*+ D'^vr^ decay from the opposite hemi- 
sphere relative to the rest of the decay products, refiecting it through the orgin, and then 
using it in the calculation of the invariant mass. 

A significant fraction of the sample of D*"*" mesons were only partially reconstructed. 
These mesons produce an enhancement in the AM spectrum very similar to the true 
signal. Only a few of the events are present in the 3-prong sample. They are more 
important in the 5-prong tagged events where a clear tail is visible in the AM distribution 
for values above 0.145 GeV (see Fig. |^d). Since such events originate from D*+ decays, 
they can still be used in the analysis. 

•^Throughout this paper charged conjugate modes are always imphcitly included. 

^In this paper any reference to the scaled energy a; of a D*+ candidate is to be understood as being 
the calculated energy of the D*+, , obtained from the reconstructed tracks, without correcting for 
missing or wrongly associated tracks, divided by the beam energy -Ebeam- 
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In all, 27 005 D*"*" candidate events were selected and the background was estimated to 
be 11 366 ± 107 events, where the error on the number of background events is statistical 
only. All of the LEP-1 data recorded by OPAL, 4.4 million events, was used for the D*"*" 
event selection. The candidate events are composed of three components: genuine D*+ 
mesons from b events, genuine D*+ mesons from c events and combinatorial background 
which is a mixture of uds, c and b quark events. Genuine D*+'s from uds quark events 
can only occur in events where a gluon splits into two heavy quarks. The possibility of 
D*+'s from this source was neglected because these events are highly suppressed due to 
the high xd*+ cut. Using Monte Carlo simulation the total contribution of this source to 
the number of D*"*" candidates was found to be (0.2 ± 0.1)%, where the error quoted is 
due to Monte Carlo statistics. 

In Ref. the fraction of genuine D*"*" mesons originating from c events was deter- 
mined by OPAL to be /° = 0.774 ± 0.023, where the error is the combined statistical 
and systematic error. The fraction of genuine D*"*" mesons from b events is given by 

/r = i-/r- 

The fractions of uds, c and b quark events in the combinatorial background, as deter- 
mined from Monte Carlo, were 0.584±0.009, 0.238±0.009 and 0.178±0.009, respectively. 
The errors are a combination of statistical errors due to finite Monte Carlo statistics and 
a systematic error accounting for the overall quality of the background estimation pro- 
cedure. Based on a Monte Carlo study |2^ the background estimate was found to be 
accurate to within 1%, and an additional 1% error was therefore assigned to the number 
of background events. Since that study was flavour blind, a conservative approach was 
taken here where the entire error was assigned in turn to each flavour component in the 
combinatorial background. 

The overall flavour composition of the c-tag event sample is presented in Table |I| along 
with the combined statistical and systematic errors. The efficiency for tagging c events 
was about 2.0%. 



6 Event shape observables 

For each flavour tagged event sample described above (uds-tag, c-tag and b-tag) the 
distributions of the event shape variables 1-Thrust (1 — T), Heavy Jet Mass scaled by 
the centre-of-mass energy (MH/y/s), Wide Jet Broadening (Bw), the ycut at which an 
event changes from being a 2-jet event to being a 3-jet event (^23) determined using 
the Durham jet finder and the C-parameter (C) were determined. The definitions for 
these observables are given in |T^,^ and the references therein. These quantities were 
measured using all tracks and electromagnetic clusters which satisfied the selection criteria 
described in Sect. ||. 

In order to extract the values of a^/a"'^^ and a^/a"'^^, the measured event shape dis- 



tributions were fitted using a QCD analytic calculation Since the QCD calculation is 
only valid for event shape distributions determined from final-state partons, hadronization 
effects caused by the transformation of final-state partons into hadrons, which are exper- 
imentally accessible, must be taken into account. This was done by applying correction 
factors to the analytic predictions. 
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The measured event shape distributions must be corrrected for experimental effects 
which distort them. These effects include finite detector resolution, initial-state photon 
radiation and biases introduced by the flavour tagging methods. In the text below de- 
scribing the correction procedure the term detector level is used to refer to distributions 
determined using the measured tracks and electromagnetic clusters and hadron level to 
refer to these distributions corrected for detector resolution, initial state radiation and bi- 
ases introduced by the flavour tagging methods. The comparisons between the measured 
distributions and the QCD predictions were performed at the hadron level. 

The event shape variable Total Jet Broadening (-Br) [0 was initially considered to be 
included in this analysis. It was found, however, that for the distribution of Bt measured 
with b quark events, the size of the hadronization corrections were greater that 20% over 
the entire distribution. Therefore this variable was dropped from further consideration. 



6.1 Correction procedure 

In the first step the selected events were corrected for distortions caused by finite detector 
resolution using an unfolding matrix. The unfolding matrix was constructed using a Monte 
Carlo data sample including initial-state radiation, full detector simulation and subjected 
to the same event selection criteria that were applied to the data. The Monte Carlo events 
which pass the selection criteria were used to calculate a correction matrix M^~*"^. The 
element M'^~^°'^{yi,yj) gives the probability that an event shape observable y measured 
at the detector level and located in bin i of its corresponding distribution, has migrated 
from bin j on the hadron level. A matrix was computed for each event shape observable 
and flavour-tagged event sample, "q-tag", where q-tag was either uds-tag, c-tag or b-tag. 

In the next step bin-by-bin corrections are used to correct the data for biases intro- 
duced by the flavour tagging methods, event acceptance and the effects of initial-state 
photon radiation. Defining G'^{yi) to be the number of events of flavour q in the un- 
tagged Monte Carlo sample and H'^~^°-^{yi) the number of events with the tag applied, 
the correction factor K'^~^°'^ {iji) for the i-th bin is given by 



The factors taken from Table |I], are the fractions of events of flavour q in each 

flavour-tagged event sample. The terms G'^{yi) are normalized to the total number of 
events of flavour q. 

Fig. D shows, as an example, the size of the bin-by-bin corrections K'^~^"'^ for the ob- 
servable By/ measured using each flavour-tagged event sample. The size of the corrections 
for Bw were typical of the other event shape observables studied. One can see for the 
distributions measured with the uds-tag and b-tag event samples that, within the chosen 
fit range, the size of the corrections are about 10%. (An explanation of how the fit range 
was chosen is presented in the next section.) In contrast, the size of the corrections for 
the c-tag event sample were an order of magnitude larger. This is due to a kinematic bias 
introduced by requiring the events in the c-tag sample to contain a high-a; D*^ meson. 

The number of events in bin i corrected to the hadron level is given by 
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j 

where N'^~^"'^{yj) is the uncorrected number of events in the j-th bin of the event shape 
distributions in question. This correction procedure does not depend on the values of 
a^^^, a'g/a^'^^ or a^/a^^^ in the Monte Carlo samples. 

The resulting distributions normalized to the total hadronic cross section are given by 

where Ayi is the width of the i-th bin and y corresponds to one of the five event shapes 
studied. 

7 Fit procedure 

The value of a^*^'^ and the ratios a^/a^*^*^ and a^/a'^'^^ were determined by simultaneously 
fitting theoretical predictions for a particular event shape observable to three hadron-level 
event shape distributions: one determined with the uds-tag event sample, one with the 
c-tag event sample and one with the b-tag event sample. The theoretical predictions were 
given by a linear combination of the theoretical predictions for each of the three tagged 
flavours: 



^ J \ q—tag,th / ^ j \ uds,t/i / ^ j \ ^^^^ 



\(^tot dy / ' \otot dy / yotot dy / 

+ fV'R{yf{—'^T' . (4) 

\(ytot dy/ 

The coefficients are the flavour fractions given in Table |I| and (l/cr^ot ■ dcr/dy)"'^^'*'^, 

{1/atot ■ da/dy)'^' and {1/atot ■ da/dy)'''*'^ are the theoretical predictions for an event 
shape observable y measured with a sample of uds, c and b quark events, respectively. 
The factors R{yY correct the theoretical prediction for a particular tagged flavour q for 
hadronization effects so that the theoretical predictions, valid for flnal state partons, can 
be compared directly with the measured distributions corrected to the hadron level. 

The hadronization correction factors R{yY were computed from JETSET 7.4 using the 
parton shower option by taking the ratio of an event shape distribution at the hadron level 
to the same distribution at the parton level. Here the parton level is deflned by the cut-off 



Qq of the QCD shower in JETSET which is set to 1.9 GeV |jT9[. After the termination 
of the parton shower the partons are transformed into hadrons using string hadroniza- 
tion |2^. For this process we choose a hybrid scheme for the longitudinal fragmentation 



function where light quarks are treated with the symmetric Lund fragmentation function 
and charm and bottom quarks according to the model of Peterson et al. [EO]. The c and 
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b quark masses in JETSET were set to their default values of 1.35 GeV and 5.0 GeV, 
respectively. 

In Equation (^), the differential cross section of a generic observable for massless 
quarks, normalized to the total hadronic cross section is given by []30| 



, \ uds,th 

1 aa\ 




The coefficients dA^'^Ydj/ and dB'^'^^'/dy are the 0{as) and QCD coefficients, re- 

spectively, atot is the one-loop cross section for the process e"'"e~ hadrons and Pq is the 
coefficient of the QCD beta function for one-loop [^. The renormalization scale n can 
be related to the e"'"e~ centre-of-mass energy by 



fi = Xf,- Ecm , (6) 

where is the renormalization scale factor. The coefficients dA^'^^/dy and dB^'^^/dy were 
obtained for each event shape observable by integrating the O(a^) matrix elements in ||30|| . 
In Equation (|^), in addition to the approximation of massless quarks, the simplifying 
assumption was made that as is the same for up, down and strange quarks. 

The terms {l/<7tot ■ da/dyY'^^ and {l/<7tot ■ da/dy)^'^^ in Equation (|) were both given 
by an Ola"^) expression for massive quarks. This expression has only recently been made 
available [|lT|. The result of this calculation was implemented in a FORTRAN program 



named ZBB4 [^] analogous to the program EVENT which was used to integrate the 
matrix elements of the massless calculation. ZBB4 was run separately for c quark 
and b quark events in order to calculate the 0{as) and 0{a'^) coefficients for the massive 
calculation. This calculation was performed in the pole mass scheme and the c and b 
quark pole masses were set to 1.35 GeV and 5.0 GeV, respectively. 

The differential cross section of a generic observable y, for massive quarks, normalized 



to the total hadronic cross section is given by ||3T 



1 da\^'*^ dAQ fafifiY 



atotdyl dy \ 2tx 



+ 27r/3olog(a;2 



^1/ 




where Q is either c or b and mg is the energy scale corresponding to the heavy quark 
pole mass. 

When performing the fit using Equation (|^) we make the substitution = a^'^^ ■ 
al/a'f'' and = ^ ■ a^/af ^ This substitution enables , al/a'f'' and a^Jaf"" 
to be determined directly as free parameters in the fit and allows correlations between 
these variables to be properly taken into account. In tests with Monte Carlo events it was 
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verified that this fitting procedure was sensitive to changes in for c and b quarks with 
respect to uds quarks. 

The fit ranges for each observable were determined by the range of y where the parton 
to hadron level corrections were below 10% and the resulting x^/d.o.f. remained small 
(~ 5 for Bw and ~ 2 — 3 for the rest). The remaining variations of the fit results due to 
the choice of the fit ranges were taken as systematic uncertainties. It was also checked that 
within the fit range the 0{a1) calculation and the JETSET parton shower model were in 
agreement with each other. The fit ranges used for this analysis are slightly smaller than 
the fit ranges used in previous OPAL studies fl^,^ which used 0{a1) calculations in 
fits to global event shapes. This is because the hadronization corrections for distributions 
of event shape variables calculated with b quark events are generally larger than the 
corrections for distributions calculated with an inclusive event sample. Consequently, the 
range of the distribution in which the corrections are below 10% is smaller. 

The fit ranges, and the results of the fits with shown in Table ^. The 

values for a^'^^ obtained appear large when compared with measurements performed 
using resummed NLO calculations [l^. This is due to the fact that fixed order QCD 
calculations, valid only to 0{a1), were used here. 

To estimate the error in our results due to the choice of the renormalization scale we 
adopted a method used in another OPAL analysis which used 0{a1) fits to global 
event shapes. We performed fits with fixed equal one, and fits in which was an 
additional free parameter. The average of the two results is taken as our main result, 
and half their difference as the "scale uncertainty". The result of fits with free 
parameter are presented in Table ^ One sees that the values of a"'^^ decrease for all event 
shapes and that the values of x^/d.o.f. show a significant improvement. This strong scale 
dependence is typical of determinations of where fixed order calculations are used. 
It was checked that the values of a^"^^ determined here were consistent with previous 
OPAL measurements of as [17, 2^. The possibility of a different scale dependence for 



heavy quark events, due to mass effects present in higher order terms, was investigated 
by introducing a separate parameter for each quark flavour. The effect was found to 
be negligible within the statistical precision of the measurement. 

In Figs. ^ through |] the results of fits with plotted along with the measured 

distributions corrected to the hadron level. The fit ranges used are indicated by the arrow 
on each figure. It can be seen that within the fit ranges the agreement between the fit 
results and the data are good. In Fig. |^ one sees a small disagreement between the fit 
results and the data for values of 1 — T > 0.25. When the fits to 1 — T were repeated 
with free parameter, the agreement between the fit results and data were good 

over the entire fit range. 



8 Systematic uncertainties 

The main result was obtained using the default selection and correction procedure de- 
scribed above. The systematic uncertainties were divided into two groups: experimental 
and theoretical uncertainties. Each uncertainty was estimated by modifying details of 
the event selection and correction procedure and repeating the analysis. The difference 
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between the results obtained with the standard analysis and the results obtained with the 
analysis corresponding to each variation were taken (unless otherwise noted) as symmetric 
systematic errors. In the case where a parameter was varied above or below its nominal 
value, the largest deviation from the main result was taken as a symmetric systematic 
error associated with that parameter. Finally, the systematic errors for each type of vari- 
ation were added in quadrature. The systematic uncertainties investigated are described 
below and the effects they had on the ratios a^/ a'^'^^ and a^/a^'^^ are presented in Tables ^ 
and |[ 



8.1 Experimental systematic uncertainties 

• The error due to the uncertainty on the flavour composition of the flavour-tagged 
event samples was evaluated by varying the flavour fractions in Table |l| within their 
errors. 

• For the central result event shape observables were measured using all charged tracks 
and electromagnetic clusters. To evaluate the relative response of the central track- 
ing and the electromagnetic calorimeter the measurements were performed again 
using electromagnetic clusters only and charged tracks only. 

• The homogeneity of the response of the detector in the endcap region was checked 
by restricting the analysis to the barrel region of the detector, requiring the thrust 
axis of accepted events to lie within the range | cos^Thl < 0.7. This only modifies 
event shape distributions measured with the c-tag event sample because the uds-tag 
and b-tag event samples are already restricted to the barrel region. 

• The minimum number of accepted charged tracks was increased from 5 to 7 in order 
to further suppress background from r^r~ events and two-photon interactions. 

• The dependency of the result on the chosen fit range was evaluated by adding or 
subtracting one bin from the lower and upper end of the fit range. Each new fit 
range was treated as a separate systematic variation. 

• In the charm event selection the requirement that X£i*+ > 0.4 leads to a large 
kinematic bias in the c-tag event sample. The uncertainty due to this aspect of the 
D*~^ selection was evaluated by varying the xci*+ cut between 0.3 and 0.5. 



8.2 Theoretical systematic uncertainties 

• The fragmentation of the final-state partons into hadrons was modelled by JETSET 
where the longitudinal fragmentation function for light quarks was treated by the 
symmetric Lund model and for heavy quarks by the model of Peterson et al ||2^ . 
The parameters for each model were determined from a fit to OPAL data on global 



event shapes [|T|]. This fit yielded a value of 6 = 0.52 ± 0.04 (PARJ(42)) for the 
Lund model and = 0.031±0.011 (PARJ(54)) and ei, = 0.0038±0.0010 (PARJ(55)) 
for the Peterson model. To evaluate the systematic error due to uncertainties in the 
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JETSET hadronization model these parameters were varied independently within 
one standard deviation of their optimised values. 



In JETSET, the parameter Qo (PARJ(82)) is the virtuality cut-off of partons. It 
determines the boundary between the pertubative QCD and hadronization phases 
and is essentially arbitrary. The optimum value to describe the OPAL data was 



determined to be Qq = (1.90 ± 0.50) GeV |]T9[ . This parameter was varied within 



one standard deviation of its optimised value. 

The width of the transverse momentum distributions of quarks and antiquarks pro- 
duced in the fragmentation process is determined by the parameter cjq = (0.40 ± 
0.03) GeV (PARJ(21)). This parameter was varied within one standard deviation of 
its optimum value. 

The renormalization scale uncertainty was estimated by performing fits with 
fixed equal one, and fits in which is an additional free parameter. The average 
of the two results is taken as our main result, and half their difference as the scale 
uncertainty. 

To evaluate the dependence of this analysis on the choice of using JETSET to cal- 
culate the hadronization corrections, the analysis was repeated using ARIADNE 
4.08 instead of JETSET. The ARIADNE parton shower is based upon a colour 



dipole model and provides an alternative to the Lund parton shower model in JET- 
SET. ARIADNE employs the same fragmentation model as JETSET for the subse- 
quent hadronization. The difference between the results using JETSET hadroniza- 
tion corrections and ARIADNE corrections were used as an estimate of this system- 
atic error. 

The error due to uncertainties in the c and b quark masses was evaluated by varying 
the quark masses in the ranges given in and repeating the integration of the 
matrix elements corresponding to massive quarks for each variation. The value for 
the b quark pole mass was varied from 4.5 to 5.5 GeV/c^; the c quark pole mass was 
varied from 1.2 to 1.9 GeV/c^. 



9 Results 

The results of the fits to determine a^/a^*^*^ and a^/a'^'^" for each of the event shape 
observables studied are presented in Tables ^ and || and summarised in Fig. |^. The values 
quoted are the average of fits with = 1 and with free parameter. The top half of 

the Fig. 1^ shows the result of determining a^/a"'^** and the lower half a^/a^'^^. The vertical 
bars on the error bars show the size of the statistical error and the full error bar is the total 
error, which is the sum of the statistical, experimental systematic and theoretical errors 
added in quadrature. Also shown is the weighted mean of the results obtained with the 
six observables. The weights were given by the reciprocal of the square of the total error 
on the ratios a^/a"*^^ and a^/a^'^^ given in Tables |^ and ^ The statistical error on the 
mean was calculated taking into account correlations between the five observables. The 
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correlation matrix was calculated from 100 Monte Carlo event samples. The systematic 
uncertainty on the weighted mean was determined from the change in the mean that 
occurred when each systematic check was applied to all the event shape distributions 
simultaneously. The mean values of a^/a"'^^ and a^/a'^'^^ were determined to be 



al/af^ = 0.997 ± 0.038 (stat.) ± 0.030 (syst.) ± 0.012 (theory) 
a^Jaf' = 0.993 ± 0.008 (stat.) ± 0.006 (syst.) ± 0.011 (theory) . 

In both cases the results are consistent with unity, indicating flavour independence of as- 
As shown in Table ^ the largest experimental systematic error on the results obtained 
from the variable ?/23 originated from moving the lower bound of the fit range from 1/23 = 
0.015 to 0.025. It was found that this relatively large variation in the result was due to a 
small discrepancy between the data and the Monte Carlo, of the order of a few percent, 
for values of ?/23 < 0.03. If one excludes the results obtained from ?/23 from the calculation 
of the weighted means, the mean value obtained for a\/a^'^^ remains unchanged within 
the experimental precision, while the mean value of a^/a^'^'^ changes from 0.997 to 0.990. 

The size of the mass effect for c and b quarks and its relevance to the measurement of 
as are also of interest. From phase space considerations one can estimate the size of the 
mass effects for 3-jet observables. The ratio of the phase space of two massive quarks and 
a gluon to the phase space for three massless particles is 1 + 8(Mq/Mz)^log(Mq/Mz) P^j. 
This represents a 7% effect for M^ = 5 GeV/c^ and 0.7% for M, = 1.35 GeV/c^. 

Fig. Ha shows the ratios al/a^'^^ and a^/a"'^^ determined using 0{a1) massless calcu- 
lations for both uds and c quarks and 0{a1) massive calculations for b quarks. Again 
we show the average of fits with = 1.0 and a;^ as a free parameter. One sees a small 
systematic shift of al/a^"^^ with respect to the fits with the massive calculation, although 
the large statistical error on the measurements for c quarks makes a definitive statement 
on the exact size of the mass effect difficult. In Fig. ||b, a^/a'^'^^ and a^/a^^^ were de- 
termined using massless calculations for uds and b quarks and massive calculations for c 
quarks. In this case a'^/a^'^^ shows a large systematic shift of order 5 to 7%. It should 
be noted that the values of a^'^^ and al/a^'^^ remained unchanged within their statistical 
errors with respect to the results presented in Fig. 0, therefore the change in a^/a^'^^ 
can be attributed entirely to changes in a^. This sensitivity to the b quark pole mass 
( Mb ) suggests that this effect could be exploited to measure Mb itself, by assuming 
flavour independence of as and fitting for Mb. Monte Carlo studies show that the value 
of Mb determined in this way depends strongly on the input value of Mb used in JETSET 
to calculate hadronization corrections, thus making a determination of Mb in this way 
problematic. 

As expected, fits to Mh do not exhibit the mass effect like the other variables do. Here 
the quark mass fixes a lower bound of the Mh distribution since the invariant mass of the 
jet cannot be less than the mass of the quark that originates the jet. For Mb = 5.0 GeV 
the lower bound lies at about Mu/y/s = 0.05, which is well outside of the fit range. 
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10 Conclusion 



We have presented a test of the flavour independence of the strong couphng constant for 
charm and bottom quarks with respect to hght (uds) quarks. This analysis was based on 
a sample of hadronic decays of the Z° resonance recorded by the OPAL detector at LEP. 
The global event shapes 7/23, 1 — T, Mh, Bw and C were used to measure as in three 
flavour tagged event samples (uds, c and b). The event shape distributions were fitted 
by 0{a^g) calculations of jet production taking into account mass effects for the c and 
b quarks. The ratios a^g/a^'^^ and a^/a^'^^ were both found to be consistent with unity, 
indicating the flavour independence of ccg. 

The measurement of a^g/a^^^ achieved a precision of 5%. The relatively large statistical 
error was due to the low efficiency for tagging c quark events (r^ 2%) which relied on 
finding D*+ mesons. The experimental error was dominated by varying the cut on the 
scaled energy of the D*+ mesons which was required to reduce the background from c 
quark events coming from cascade decays of b hadrons. The measurement of a^/a"*^^ 
achieved a 1.5% precision, the error being dominated by the theoretical systematic error. 
The largest theoretical errors were due to uncertainties in the b quark mass and the 
renormalization scale factor . 

In addition, we have presented a study of the effect of heavy quark masses on global 
event shape variables. It was observed that the values of as determined from 1/23, 1 — T, 
Bw and C were reduced by 5 to 7% when these event shapes were measured with a sample 
of b quark events and a massless QCD calculation was used. The shape variable Mh was 
found to be insensitive to the b quark mass. 
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Tables 





uds fraction 


charm fraction 


bottom fraction 


uds-tag 


0.862 ±0.010 


0.119 ±0.003 


0.019 ±0.011 


c-tag 


0.246 ± 0.005 


0.548 ±0.018 


0.206 ±0.017 


b-tag 


0.008 ± 0.002 


0.028 ± 0.004 


0.964 ±0.019 



Table 1: The fraction of uds, c and b quark events in each flavour-tagged event sample. 
The errors are the combined statistical and systematic errors. 



cut 


3 prong 


scmilcptonic 


satellite 


5 prong 




0.4-1.0 


0.4-1.0 


0.4-1.0 


0.5-1.0 


Mdo [GeV] 


1.79-1.94 


1.20-1.80 


1.41-1.77 


1.79-1.94 


AM [GeV] 


0.142-0.149 


0.140-0.162 


0.141-0.151 


0.142-0.149 



Table 2: List of cuts used in the D*+ reconstruction. Note that both the scaled energy 
xd*+ and the mass difference AM are effective quantities, calculated from the recon- 
structed tracks only. The exact meaning of the different quantities is explained in the 
text. 
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Event shape 


Fit Range 


^uds 






xVd.o.f. 


1/23 


0.015 - 0.185 


0.1263 ±0.0006 


1.017 ±0.034 


1.004 ±0.007 


1.4 


1 -T 


0.11 - 0.29 


0.1451 ±0.0011 


0.967 ±0.055 


0.993 ±0.011 


2.2 


Mh 


0.26 - 0.42 


0.1344 ±0.0009 


1.004 ±0.047 


1.008 ±0.009 


2.5 


Bw 


0.09 - 0.21 


0.1313 ±0.0007 


1.006 ±0.039 


1.003 ±0.008 


4.4 


C 


0.35 - 0.67 


0.1424 ±0.0010 


0.986 ±0.050 


0.982 ±0.011 


1.6 



Table 3: Results of the fits with the renormalization scale factor fixed: = 1. The 
errors are statistical. 



Event shape 








•^^ 


xVd.o.f. 


y23 


0.1188 ±0.0010 


1.008 ±0.030 


0.991 ±0.006 


0.266 ±0.066 


1.0 


1-T 


0.1163 ±0.0017 


0.963 ±0.036 


0.979 ±0.007 


0.083 ±0.017 


0.5 


Mh 


0.1195 ±0.0006 


0.990 ±0.037 


0.992 ±0.008 


0.088 ±0.021 


1.2 


Bw 


0.1191 ±0.0006 


0.995 ±0.033 


0.995 ±0.007 


0.086 ±0.011 


1.8 


C 


0.1125 ±0.0014 


0.973 ±0.032 


0.984 ±0.010 


0.041 ±0.014 


1.1 



Table 4: Results of the fits with the renormalization scale factor free parameter. 

The errors are statistical. The fit ranges are given in Table |^. 
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Event shape 


2/23 


1 -T 




Bw 


C 


Mean 




1.013 


0.965 


0.997 


1.001 


0.979 


0.997 


Statistical error 


±0.032 


±0.045 


±0.042 


±0.036 


±0.041 


±0.038 


Flavour composition 


±0.002 


±0.002 


±0.001 


±0.001 


±0.002 


±0.001 


Clusters only 


+0.015 


+0.011 


+0.014 


+0.014 


+0.033 


+0.017 


Tracks only 


+0.009 


+0.012 


+0.037 


+0.006 


±0.001 


+0.012 


COS^Th < 0.7 


+0.005 


-0.008 


+0.027 


+0.010 


+0.022 


+0.011 


Nch > 7 


< 0.001 


+0.001 


+0.001 


+0.001 


+0.001 


+0.001 


Fit range 


±0.008 


±0.016 


±0.010 


±0.013 


±0.003 


±0.004 


xd*+ range 


±0.009 


±0.045 


±0.020 


±0.008 


±0.037 


±0.019 


Experimental syst. 


±0.022 


±0.051 


±0.053 


±0.024 


±0.055 


±0.030 


Lund frag, b 


±0.003 


±0.001 


±0.002 


±0.002 


±0.001 


±0.001 


Peterson frag. Sc 


±0.010 


±0.003 


±0.002 


±0.001 


±0.003 


±0.003 


Peterson frag, 


±0.005 


±0.004 


±0.004 


±0.004 


±0.003 


±0.004 


Qo 


±0.006 


±0.001 


< 0.001 


±0.003 


±0.003 


±0.003 




±0.004 


±0.003 


±0.014 


< 0.001 


±0.005 


±0.002 


Scale uncertainty 


±0.004 


±0.002 


±0.007 


±0.005 


±0.006 


±0.007 


ARIADNE 


-0.008 


< 0.001 


-0.001 


-0.004 


-0.001 


-0.004 


c quark mass 


±0.008 


±0.005 


< 0.001 


±0.005 


±0.006 


±0.005 


b quark mass 


±0.001 


< 0.001 


< 0.001 


±0.001 


< 0.001 


< 0.001 


Theoretical error 


±0.018 


±0.008 


±0.016 


±0.010 


±0.011 


±0.012 


Total error 


± 0.043 


±0.069 


±0.070 


±0.044 


±0.069 


±0.050 



Table 5: Systematic errors on the ratio a'^/a^'^''^. Where a signed value is quoted, this 
indicates the direction in which the ratio changed with respect to the default analysis 
when a certain feature of the analysis was changed. See the text for a description of the 
systematic errors. 
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Event shape 


2/23 


1-T 






c 


Mean 




0.997 


0.986 


1.000 


0.999 


0.983 


0.993 


Statistical error 


±0.007 


±0.009 


±0.008 


±0.008 


±0.010 


±0.008 


Flavour composition 


±0.001 


±0.001 


< 0.001 


< 0.001 


±0.001 


< 0.001 


Clusters only 


-0.004 


±0.001 


+0.001 


+0.001 


±0.001 


< 0.001 


Tracks only 


-0.006 


+0.002 


+0.003 


+0.001 


-0.001 


-0.001 


COS^Th < 0.7 


< 0.001 


-0.001 


+0.002 


+0.001 


+0.002 


+0.001 


Nch > 7 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


Fit range 


±0.016 


±0.008 


±0.009 


±0.003 


±0.003 


±0.005 


xd*+ range 


±0.001 


±0.004 


±0.002 


±0.002 


±0.003 


±0.002 


Experimental syst. 


±0.017 


±0.010 


±0.011 


±0.004 


±0.005 


±0.006 


Lund frag, b 


±0.001 


±0.005 


±0.003 


±0.001 


±0.005 


±0.001 


Peterson frag. Cc 


±0.003 


±0.007 


±0.005 


±0.002 


±0.001 


±0.003 


Peterson frag, ej, 


±0.004 


±0.004 


±0.003 


±0.001 


±0.004 


±0.002 


Qo 


±0.004 


±0.008 


±0.003 


±0.003 


±0.005 


±0.003 




±0.003 


±0.009 


±0.005 


< 0.001 


±0.001 


±0.001 


Scale uncertainty 


±0.006 


±0.007 


±0.008 


±0.004 


±0.001 


±0.005 


ARIADNE 


+0.004 


-0.001 


-0.001 


+0.008 


-0.002 


+0.002 


c quark mass 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


< 0.001 


b quark mass 


±0.012 


±0.008 


±0.005 


±0.010 


±0.010 


±0.008 


Theoretical error 


±0.016 


±0.018 


±0.013 


±0.014 


±0.013 


±0.011 


Total error 


±0.025 


±0.022 


±0.019 


±0.016 


±0.017 


±0.014 



Table 6: Systematic errors on the ratio a^/a"'^''. Where a signed value is quoted, this 
indicates the direction in which the ratio changed with respect to the default analysis 
when a certain feature of the analysis was changed. See the text for a description of the 
systematic errors. 
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Figures 




2 4 6 8 10 



Figure 1: Distribution of number of tracks per event (N^jg) with b/ai, > 2.5. The 
points are data and the sohd histogram is Monte Carlo. The Monte Carlo is also shown 
decomposed into contributions from uds, c and b quark events. 
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Figure 2: Distribution of the mass difference AM for four different D*"*" decay cliannels: 
(a) 3-prong decay mode, (b) the two semileptonic modes combined, (c) the sateUite decay 
mode, and (d) the 5-prong decay mode. The arrows indicate the selected signal regions. 
The points with error bars are the signal candidates and the solid histograms are the 
background estimator distributions. 
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Figure 3: Bin-by-bin correction factors applied to the distribution of B^r for (a) the uds- 
tag event sample, (b) the c-tag event sample and (c) the b-tag event sample. The fit 
range for this observable (see Sect. |^ is indicated by the arrow on each figure. 
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Figure 4: The points with error bars are the distributions of 1/23 and 1 — T measured in 
three flavour-tagged event samples (uds-tag, c-tag and b-tag) corrected to the hadron level. 
The curves show the result of simultaneous fits of 0{a1) calculations to the distributions 
measured in each flavour-tagged event sample. The renormalization scale a;^ was flxed 
equal to one. The flt range for each distribution is indicated by the arrows. 
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Figure 5: The points with error bars are the distributions of Mh and measured in 
three flavour-tagged event samples (uds-tag, c-tag and b-tag) corrected to the hadron level. 
The curves show the result of simultaneous fits of 0(a;^) calculations to the distributions 
measured in each flavour-tagged event sample. The renormalization scale was flxed 
equal to one. The flt range for each distribution is indicated by the arrows. 
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Figure 6: The points with error bars are the distribution of the C parameter measured in 
three flavour-tagged event samples (uds-tag, c-tag and b-tag) corrected to the hadron level. 
The curves show the result of simultaneous fits of 0{al) calculations to the distributions 
measured in each flavour-tagged event sample. The renormalization scale Xfj, was fixed 
equal to one. The fit range for each distribution is indicated by the arrows. 
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Figure 7: The ratios a'^/a'^'^^ and a^/a^'^^ determined for each of the event shapes studied. 
The weighted mean of the results derived for the five event shape observables is also shown. 
The inner error bar represents the statistical error while the full error bar represents the 
statistical error and the systematic error added in quadrature. 
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Figure 8: The ratios a'g/a^'^^ and a^/a"'^^ determined in two different ways: (a) where 
the uds and c quarks were treated as if massless and the b quark was considered to be 
massive and (b) where the uds and b quark were treated as massless and the c quark was 
considered massive. The errors shown are statistical only. 
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